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ABSTRACT 


A one-dimensional mixed layer model (Garwood, 1976, 1977) is used to 
Simulate the mixed layer depth and temperature observed during POLEX, a 
component of the North Pacific Experiment. POLEX occurred during Janu- 
ary-February 1974 with mid-ocean observations (35°N, 155°W) of tempera- 
ture and salinity made from R/P FLIP, while FLIP was under free-drift 
conditions. The results of this simulation show that apparent horizon- 
tal advection, due to the drift of R/P FLIP, was important in the heat 
and salinity budgets, but that vertical mixing and mixed layer depth 
changes were controlled primarily by the one-dimensional response of the 
turbulent kinetic energy budget to the local atmospheric forcing. Oc- 
casionally, surface salinity flux, due to large precipitation rates, can 
Significantly alter (on the order of tens of meters) the thickness of 
the ocean surface turbulent boundary layer, with a demonstrated decrease 


in thickness during a single period of strong precipitation. 
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NOTATION 


Buoyancy 

Turbulent buoyancy flux (cm’sec 3) 
Degrees Centigrade 

Non-dimensional drag coefficient 

Heat capacity of sea water (cal gm7!c7!) 


~e Sect!) 


Rate of evaporation (gm cm 
Vapor pressure of air 
Saturated vapor pressure of marine air 
Coriolis parameter 

Acceleration due to gravity 

Turbulent boundary layer depth (modeled) 
Mixed layer depth (modeled) 

Observed mixed layer depth 

Meters 

Rate of precipitation (gm cm™“sec™ |) 
Upward surface heat flux (cal cm *sec™!) 


-2 


Back Radiation cal cm sec’ |) 


( 
Latent heat flux (cal cm “esac! ) 
Sensible heat flux (cal cm ~e5e¢ ) 
Clear sky radiation (cal a ) 
Net surface heat flux (cal eTacee ) 
Salinity (0/00) 
Turbulent salinity flux (0/00 cm secre!) 


Temperature 











Temperature in the air 

Sea-surface temperature 

Turbulent temperature flux (mwatts cm7*) 
Time (hours) 


Frictional velocity at ocean surface 


Mean wind speed (m sec’ !) 


Vertical turbulent velocity flux 
Velocity vector (cm sec’) 

Velocity vector due to ocean currents 
Velocity vector due to Ekman flow 
Velocity of current due to quasi-geostropic flow 
Velocity due to the drift of R/P FLIP 
Vertical component of turbulent velocity 
Depth (meters) 

Depth increment 

Expansion coefficient for heat 

Expansion coefficient for salt 

Density of sea water (gm cm™>) 


Density of air (gm m7 3) 


Surface stress (dynes m2) 
Temperature 
Represent the total instantaneous value 


Denotes the vertical average 


Increment of change in a variable 
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I. INTRODUCTION 


A. PURPOSE OF STUDY 

The objective of this research was to investigate the conclusions 
which were reached by Simpson and Paulson [1979] from ocean data ob- 
tained as part of the Pole Experiment (POLEX), a component of the North 
Pacific Experiment (NORPAX). POLEX was part of the first process ori- 
ented NORPAX experiment. The experiment was named POLE to indicate that 
the horizontal extent of the sampling was small compared to the largest 
scales investigated in NORPAX. Simpson and Paulson concluded: '"Qne- 
dimensional mixed-layer deepening models failed to predict the mixed- 
layer depth and temperatures observed during POLE. Horizontal advec- 
tion, as evidenced from the salinity maximum frequently occurring at the 
bottom of the mixed layer and other near-surface changes in salinity and 
temperature not associated with local surface forcing, are responsible 
for the failure". 

This study will investigate the one-dimensionality of the POLE 
temperature and salinity data and determine if the mixed layer depth 
changes in response to the local atmospheric conditions (one-dimensiona] ) 
or if these changes are caused by two and/or three dimensional advective 
processes. To accomplish this task, a non-stationary, one-dimensiona | 


bulk model of the mixed layer originally proposed by Garwood [ 1976, 1977] 


Simpson, J. J. and Paulson, C.A., 1979: Observation of Upper Ocean 
Temperature and Salinity Structure During POLE Experiment, J. Phys. 


Oceanogr., 9, page 869. 





is expanded by including a salinity budget. A comparison of the re- 


sults of a strictly one-dimensional model of the upper ocean with obser- 
vations can be used to make an assessment of the one-dimensionality of 


the different variables which are involved in the upper ocean boundary 


layer response to the local atmospheric forcing. 

One hypothesis to be examined is that horizontal advection may be 
important for the heat and salt budgets, and at the same time vertical 
mixing and mixed layer depth may be controlled primarily by one- 
dimensional or local forcing alone. One objective of this study will 


be to determine how responsive are the variables of heat and salt in 
the upper ocean boundary layer to one-dimensional effects. A second 
objective of this study is to use the data obtained as part of POLEX to 
see 1f the one-dimensional mixed layer model proposed by Garwood will 
predict the observed mixed layer depth, even if advection dominates the 
local changes in heat and salt content in the upper ocean. 

The data used in this report contain observations made from R/P FLIP 
while under free-drift conditions. As a result of the free-drift condi- 
tion of FLIP during the data collection period, advection as seen by R/P 
FLIP will be a combination of the actual advection in an Eulerian coor- 
dinate system for the area and advection due to the drift of FLIP. 
Henceforth, the term advection will be used to mean a combination of the 
actual advection term and apparent horizontal advection due to the drift 
erent IP. 

The data were collected during the period 28 January 74 through 
14 February 74. During that time, FLIP occupied a station approxi- 
mately 800 miles north of the Hawaiian Island Chain under free drift 


conditions. The position of FLIP ranged from 35°39'N, 155 °05'W to 


10 





34°36'N, 155 25'W. The experimental site is hydrodynamically complex, 
as shown in Figure 1. The Subtropical Front is known to meander between 
32° and 35°N (Roden, 1974; Barnett,1976). The region of the trade winds 
northeast of Hawaii has air-sea fluxes of latent heat in excess of 35 
cal cm” “day”! (Wyrtki, 1965). The Subtropical water mass formed in this 
region contrasts markedly with the less saline Eastern North Pacific 
Central Water characteristically encountered north of 35°N. The Horse 
Latitudes are located only 3° of latitude to the south of the observa- 
tion area and the North Pacific Current is expected to affect the gen- 


eral hydrography of the region. 


B. MILITARY IMPORTANCE 

Mixed layer depth is an important factor when dealing with sound 
propagation in the ocean. It is used to determine such sonar detection 
variables as convergence zone distances, source level depth, and cutoff 
frequency for surface duct propagation. A model which could predict the 
change in mixed layer depth in response to the local atmospheric forcing 
could be incorporated in the existing acoustic models, for example the 
Fast Asymptotic Coherent Transmission (FACT) and the Naval Interim Sur- 
face Ship Model (NISSIM), presently being used at Fleet Numerical Ocean- 
ography Center (FNOC). A one dimensional mixed layer model would help 
the naval tactician to directly correlate the variability of the ocean 
with hisweapon system and determine the most efficient deployment of 


the naval Anti-Submarine Warfare (ASW) assets. 


C. LITERATURE REVIEW 


The vertical fluxes of heat, salt, and momentum across the sea-air 


interface are the sources of almost all oceanic motions. In the fully 
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turbulent oceanic mixed layer bounded by the sea-air interface above and 
by the dynamically stable water mass below, the vertical fluxes are 
large. Below the mixed layer, they are usually negligibly small so that 
one can decouple the mixed layer from the underlying stable, quiescent 
water mass. This homogeneity is the root of the term "slab", which is 
often used to describe the layer as depicted in Figures 2 and 3. 

When dealing with sound propagation in the ocean, the term mixed 
layer depth is sometimes taken to be that area with a uniform tempera- 
ture structure, but this is not always the case. In Figure 4 the 
bottom of the mixed layer depth (Region I) appears to be at about 80 
meters, and below there is a strong temperature gradient. However, ex- 
amination of the on profile shows that the mixed layer depth is about 30 
meters. This shallow halocline was probably caused by an earlier rain- 
fall which decreased the near-surface salinity. The assumption of ver- 
tical homogeneity in a bulk model, because of the problem associated 
with solving for the interior fluxes of buoyancy and momentum, reduces 
to the need to Know only the surface and entrainment fluxes. However, 
only small vertical gradients in these mean variables may be associated 
with large turbulent fluxes. Therefore, the slab assumption should not 
be as automatically applied to the turbulent kinetic energy budget. 

The immediate local reaction of this mixed layer to those fluxes re- 
sults in a homogeneous water column, i.e., vertical uniformity in the 
mean velocity and density. This adjustment to the density and velocity 
Structure of the surface layers of the ocean to variable fluxes has-been 
the subject of a large number of studies since Ekman's [1905] treatise 


where he originated the concept of a depth of frictional resistance for 





the upper section of a wind stressed ocean. He further suggested that 
the wind-driven current will have a similar depth, the Ekman depth. 

Much of the one-dimensional theory for the ocean surface turbulent 
boundary layer or mixed layer seems to be dependent upon the validity 
of two basic hypotheses. The first of these is that vertical mixing 
within the turbulent boundary layer and entrainment mixing at its base 
occur in response to local atmospheric forcing, i.e., the surface_wind 
stress and net heat radiated or fluxed across the sea-air interface. The 
second hypothesis is that the mechanical energy budget is the key to the 
understanding and prediction of the mixed layer dynamics. 

Most observations show that local atmospheric forcing such as an 
increase in wind speed increases the mixed layer depth and low wind speed 
and high solar radiation decreases mixed layer depth. This is less true 
with increased depth below the mixed layer. Observations by Pollard and 
Tarbel [1975] showed that horizontal coherence between the fluctuating 
wind-driven currents near the surface was large and consistent with the 
local amplitude and phase differences of the wind history. With in- 
creased depth, the horizontal coherence decreased monotonically as did 
the vertical coherence of the deeper current measurements with those 
near the surface. Kroll [1975] pointed out that the time scale of these 
Strongly inertial currents is on the order of a couple of days, depend- 
ing upon the stability of the underlying water column as well as the 
difference between the driving frequency and the local inertial fre- 
quency. The turbulent motions within the mixed layer are expected to be 
even more responsive to the local atmospheric conditions because of the 


small dissipation -time scale for the turbulent kinetic energy--minutes 


ie 





to hours, which depends upon the layer depth and turbulence intensity. 
Both the mean and turbulent kinetic energies in the mixed layer are ex- 
pected to reflect the local boundary conditions for wind stress and sur- 
face buoyancy flux. Measurements of turbulent velocity fluctuations, 
which were made by Powell et al [1975] and Jones and Kenney [1977], show 
that the root-mean-square turbulent velocity does scale with the water 
surface friction velocity, u,. Camp and Elsberry [1978 | and Elsberry 
and Camp [1978] demonstrated the role of strong synoptic scale wind 
events for deepening. Garwood and Halpern [1976] and Elsberry and Raney 
[1978] indicated the importance of diurnal heating, especially at times 
between storms. Hindcasting simulations by Thompson [1976], Mellor and 
Durbin [1975] and Camp [1976] support the importance of surface boundary 
fluxes. 

Kraus and Turner fi67] improved the one-dimensional model originally 
developed by Kraus and Rooth figei] . This was accomplished by consider- 
ing the turbulent kinetic energy budget utilizing the heat equation and 
a mechanical energy equation. Because the frictional generation of heat 
is negligible, the vertically integrated heat equation provides a rela- 
tionship for the conservation of potential energy. They parameterized 
the mechanical production rate in terms of the friction velocity, but 
neglected the viscous rate of dissipation and the effects of entrainment 
Shear production. 

The application of mechanical energy principles to the problem of 
mixed layer dynamics has caused some disagreement. An important contro- 
versy in mixed layer modeling has been centered around whether mean kin- 


etic energy or turbulent kinetic energy leads to the downward entrainment 
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buoyancy flux at the base of the turbulent boundary layer. The mean 
kinetic energy advocates, including Thompson [1976], Pollard [i977] and 
Price [1977] » argue that this mean flow energy is directly converted to 
potential energy by a mean flow instability that causes the bulk 
Richardson number to be less than some critical constant of the order 
one. Turbulent kinetic energy advocates do not disagree with the pos- 
sibility of this mechanism, but they argue it is not the dominant mech- 
anism because turbulence-generated instabilities will usually suffi- 
ciently erode the interface to preclude the occurrence of the mean flow 
instability. De Szoeke and Rhines [1975] and Yun [1978] demonstrated 
that wind-shear production of turbulent energy is ultimately the most 
important source of energy for deepening the mixed layer, if the surface 
buoyancy flux is neglected. Even when the mean flow energy is important, 
they find that the rate of conversion of this energy to potential energy 
via entrainment shear production is controlled by the constraints of the 
turbulent kinetic energy budget. The Kraus-Turner turbulent energy model 
has not proven to be entirely satisfactory, and a number of changes have 
been suggested. Stevenson 11979] examined a number of alternative para- 
meterizations for dissipation. Zeman and Tennekes [1975], Elsberry et al 
fi976], Alexander and Kim [1976] and kim [1976] added depth-dependent 
dissipation terms that helped to reduce unrealistic deepening: 

Garwood l1976, 1977 developed an ocean mixed-layer model using the 
Navier-Stokes equation of motion with the geostrophic component eliminated, 
the continuity equation in incompressible water, the heat equation from 
the first law of thermodynamics, the conservation of salt equation, and 


a linearized equation of state. The entrainment hypothesis depends upon 


15 





the relative distribution of turbulent energy between horizontal and 
vertical components and is offered as a plausible mechanism for govern- 
ing both entrainment and layer retreat. 

Garwood suggested that planetary rotation influences the dissipation 
of turbulence for deeper mixed layers and enables a cyclical steady state 
to exist on an annual basis. Furthermore, the rate of entrainment for the 
stable regime can not be a simple linear extrapolation of the unstable 
Situations. Unlike the atomospheric boundary case, most of the solar 
radiation does not penetrate the layer. Therefore, downward turbulent 
heat flux in the oceanic boundary layer is as important as the upward 
flux during the course of both diurnal and annual cycles. The non- 
linearity of the interface entrainment tendency parameter, which is 
greatest for stable surface boundary conditions, results in a modulation 
of the long-term trend of the mixed-layer depth by the diurnal component 
of surface heat flux. In this model, buoyant production is somewhat 
more efficient than shear production as a source of energy for vertical 
mixing because of its unique effect on the vertical component of the 
turbulent velocity. The buoyancy equation is generated from the heat 


and salt equations together with an equation of state, 


B= oy [l-a@-) +8 %-s,)], (1) 


and the definition for buoyancy, 
a AV 
b = g (p49 - £)/p, (2) 
AS) 20 
In equations 1 and 2 6 is temperature, S salinity and p density 


while a and B are the expansion coefficients for heat and salt, respec- 


tively, and g is gravity. The tilde represents the total instantaneous 
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value and the subscript zero denotes a representative but arbitrary con- 
stant value. The generalization of using b rather than 8 will cast the 
model equations in a form equally applicable to those situations where 
evaporation and precipitation contribute significantly to the surface 
buoyancy flux and the structure of the evolving pycnocline. The buoyancy 
equation also has a more obvious and direct role in the mechanical energy 
budget as shown in Figure 5. 

im tter i976] included salinity in the mixed-layer model of Kraus and 
Turner [1967] . He concluded that salinity is important in determining 
the density structure when applying a mixed-layer model to the entire 
ocean, but the relative importance of salinity effect on the short term 
evolution of the density profile may not be significant. Paulus 1978] 
included salinity in the mixed layer model of Camp [1976} . He concluded 
that the effects of salinity structure were noticeable only in deepening 


regimes. 


ile 





IT. DATA ANALYSIS 


A. METHODS 

The upper ocean temperature and salinity data and the meterological 
data during the NORPAX POLE experiment were obtained from Simpson and 
Paulson fi977]. The parameters required for surface boundary condition 
computations, which are used in Garwood's model, include wind speed, wind 
direction, cloud cover, sea surface temperature, air temperature (dry 
bulb), dew point and rain code. These parameters were recorded at ir- 
regular time intervals. A linear interpolation technique was used to 
obtain values for these input parameters at regular time intervals. The 
resulting interpolated data were compared with the original data. Fig- 
ures 6 through 12 show that no noticeable error was introduced into the 
parameters by the use of a linear interpolation. 

Temperature and salinity profiles as in Figures 3 and 4 were digi- 
tized by using a Houston Instrument Hi Pad Digitizer and a Tektronix 
4052 computer. The digitized profiles of temperature and salinity were 
linearly interpolated into hourly intervals. This provided for both 
initial conditions in the one-dimensional model and for fields to be 
used to compute observed mixed layer depth, h(t), and to compare with 


subsequent model computations of T(z,t) and S(z,t). 


B. SURFACE FLUXES 
The mixed layer model uses the total surface heat flux (Q,(0)); in- 
cident solar radiation, wind speed at a height d, rate of evaporation 


(E), and rate of precipitation (P) to calculate surface fluxes of 
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buoyancy (heat and salt) and momentum. The friction velocity, in air, 


is calculated using the following formulas: 


> p4Cgilg (3) 
Us = ner (4) 


where OL is the mean wind speet (m/sec), 


Cy is the non-dimensional drag coefficient (1.4x1073), 
P, is the density of air (1.25x1072gm/cm?) , and 


as the surface stress (dynes/cm°). 
The turbulent fluxes of latent heat (Q,) and sensible heat (Q,) are 
estimated using the following bulk aerodynamic formulas: 


Q, = Cy(-98E, - E,)u, (5) 


Renee > t)U (6) 
The net back radiation (Q.) is estimated from the following empirical 
formula reported by Husby and Seckel [1975]: 

Q, = 1.14 x 10°7(273.16 + T.)° (.39 - .05E,/*) (1. - .6c%) (7) 
where E. is the saturated vapor pressure of the marine air directly in 
contact with the sea surface (.98 corrects for salt defects), 

T, is the air temperature (degrees Centigrade), 


a 

E 1s the vapor pressure of air at approximately 10 m based on dew 
point temperature, 

ee is the sea-surface temperature (degrees Centigrade), and 

C is the fractional cloud cover. 


The upward heat flux is then: 
q,=0, +0, + 0, (8) 





The solar insolation, Q.> was estimated by: 
Q. = (1. = aa?) (1. - .66C°)Q, (9) 
The constants a and b are adapted from Tabata [1964] and the cubic 
cloud cover correction from Laevast [i960] . The coefficient C is the 
fractional cloud cover and ~ is the mid-day elevation angle of the sun. 
The clear sky radiation, Qo» is given by 
Qp=AptA; cosotB, sindtA,cos26+B,s in2o (10) 
from the formula developed by Seckel and Beauday [1973]. The coeffi- 
cients (Ap» etc.) were calculated by harmonic representation of the values 
presented in the Smithsonian Meteorological Tables (List, 1958) with 
o=(27/365)(+-21) where +t is the julian day of the year. Turbulent 


temperature flux, 





Tw’ (9) =-Qa_ al) 


Cc 
ap 


turbulent velocity flux, 





turbulent salinity flux, 





me (0) ~ ‘F-£)5(0) as 


and turbulent buoyancy flux, 





where a is the coefficient of thermal expansion, 
B is the coefficient of contraction of sea water due to salinity 
variations, and 


5(9) is the surface salinity value. 
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Figures 13, 14, 15 and 16 are graphs which show the hourly values of 
wind stress, tT, total heat flux at the surface, Q,, salinity filiee S?w" , 
and buoyancy flux, b'w', respectively, as calculated by the one-dimensional 
model. These figures can be compared with Figures 17 and 18 which show the 
interpolated hourly values of wind stress, t, total heat flux at the sur- 
face, Qe» buoyancy flux, b'w', and various other components which were re- 
ported by Simpson and Paulson [1977]. Wmeyvalues ef these fluxes cal- 
culated here are nearly identical to the values reported by Simpson and 
Paulson. In Simpson and Paulson's graph on buoyancy flux (Figure 17), a 
positive buoyancy flux results when ae > 0. In the one-dimensional 
model (Figure 16) a positive buoyancy flux results when Oe eiOand. — ip. 

1. Surface Wind Stress 

The wind stress calculated using the bulk approximation with a 
drag coefficient of 1.4x107> has a mean value of about 0.66 dynes cm~? 
for the experimental period. Figure 13 shows the output for the surface 
wind stress as calculated by the one-dimensional model. Analysis of 
Figure 13 shows that for periods up to as long as 4 days the wind stress 
1s considerably less than 1 dyne cm~?, and it 1s evident that only two 
periods of relatively high winds were encountered during the experiment. 
One period of high wind near the middle of the experiment (time = 200-260 
hours) with a maximum wind stress of 2.5 dynes em™?, and the other at the 
end of the experiment with a maximum wind stress of 3.3 dynes om™*. 

2. Surface Heat Flux 

Comparison of the total heat flux with net all-wave flux (Figures 
17 and 19) shows that the net all-wave flux is generally the dominant 


term in the surface heat balance. However, a two-day period of enhanced 
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latent heat flux, due largely to a significant increase in wind speed 
midway through the experiment (time = 250 hours), is evident when com- 
paring Figures 13 and 14. 

Simpson's [1977] calculation of the total heat budget for the 
experiment suggested that the ocean gained 201 cal cm™* for the time per- 
iod from 3 Feb 74 - 12 Feb 74. Simpson's heat budgets calculations indi- 
cated that ocean gained heat for 6 days, was in near thermal equilibrium 
with the atmosphere for 2 days and lost heat to the atmosphere for 2-day 
period. These results, coupled with the low wind stress, imply that the 
upper dynamics might have been dominated by net surface heating for a 
Significant part of the experiment. 

eee oalinity Flux 

Figure 15 shows the salinity term as computed by the one- 
dimensional model (equation 13) without including the precipitation term 
in the computation. The precipitation term was omitted in the calculations 
because no reasonable correlation could be made between a rain code given 
in the experiment data report (Simpson and Paulson 1977) and an amount of 
precipitation. The results of model runs which included precipitation 
(estimated by the rain code) were compared to model runs which neglected 
precipitation altogether. These model comparison runs showed no signifi- 
cant difference except during one period when precipitation was much 
larger than evaporation. Simpson and Paulson [1979 observed that over- 
all evaporation and precipitation were in near equilibrium, with precip- 


itation exceeding evaporation by 23 mg em7°. 


4. Buoyancy Flux 


Stratification in the ocean exerts a strong dynamical influence 


Ze 





on mixing since a stably stratified fluid requires work against buoyancy 
forces if mixing is to occur. For a stably stratified fluid, the density 
increases with depth and the buoyancy flux is directed downward, while 
the converse is true for an unstable, stratified fluid. In the one- 
dimensional model, surface buoyancy flux, B'w' (gy is determined by the 
air-sea transfers (equation 14). Figure 16 shows the model computed 
buoyancy flux with ORs < 0 and E > P resulting in a positive buoyancy 
flux. The net heat flux, Figure 14, dominates the surface values of 
buoyancy flux during the experiment. During daytime periods, negative 
buoyancy flux results from the domination of the solar heating at the 
Surface, while during the night periods a positive buoyancy flux is 


produced by the combined effect of net long-wave radiation and the up- 


ward turbulent fluxes of heat and moisture. 


C. ADVECTION 

The determination of the importance of advection to the heat and 
Salt budget is necessary before an assessment can be made on how one- 
dimensional is the response of the variables of heat and salt are in 
the upper boundary layer. The velocity vector V of the advection term 
(V' (VT,VS)), is separated into two components in order to determine 
which component dominates the advection term. 

1. Heat and Salinity Budgets 

Neglecting internal sources and sinks, the equation for conser- 


vation of heat or salt is: 


0 

f |aciseto0rs) + 9(7,5)'w' + a(T,S)'u' + Tar | dz=0 (Ts) 
dt a. Z 0 xX dy 

-d 


(a) (b) Ke) (d) (e) 
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where T 1S heat, S is salinity, t is time, d is depth and V is the velo- 

city vector in a coordinate system moving with the drifting platform. 

The meanings of these terms are: 

(a) Heat or salinity change with time, 

(b) Heat or salinity change due to advection in the moving coordinate 
system, assuming that when averaged over a day such features such 
as internal waves would not contribute to advection, 

(c) Vertical heat or salt flux divergence, and 

(d) and (e) Horizontal heat or salinity flux divergence. 

Terms (d) and (e) are assumed to be small because they can not be 
evaluated. Therefore, equation (15) can be rewritten in terms of the 


variable temperature: 


0 
| 2 tor | dz =- (T'w' (gy) + Tw (d) (16) 
= (f) 





Term (f) is equal to the surface heat flux Qr/PC, » where Qr 1s the net 
surface heat (Q7=Q. -Q-Q,-Q, +0.) 5 cal cm” sec’! 0 is the density of sea 
water, | gm cm73, and Cy is the specific heat of sea water, 1 cal gm co! 
Term (g) is assumed to be negligible if d is well below the mixed layer, 
ieee, d > aay The vertical integrated mean advection of temperature is 


therefore, 
0 


0 0 t 
- Je [free], [free]. | fore on 
_ Rome, tt 
d ite O pc, 


24 





The same procedure used to derive (17) may be followed to find the 
average advective change in salinity. The surface salinity flux is 
given in equation (13). Thus the vertically integrated mean advection 


of salinity is: 





0 0 0 “ 
| 7 FsFas - fi : fee + [sere (18) 
- - 
-d -d ‘d : 
tg 
At 


The magnitude of actual advection (in an Eulerian coordinate system) 
can be checked by estimating V and VT or VS. The term V'(VT, VS), is the 
advection quantity averaged over a time interval, At. The velocity vector 
V can be separated into two components, 

V=V.+V (19) 
where V_ 1S the velocity vector due to the ocean currents and V 1s the 
velocity vector due to the drift of R/P FLIP. The velocity vector due to 


the ocean currents 1s itself composed of two components: 


Me = Vg a Ve 020) 
where V_ is the velocity of the current due to the quasi-geostrophic flow 


g 
and is assumed to be parallel to the surface temperature contour. There- 


fore V:V7 1S assumed to be negligible. The order of magnitude of V., 
the velocity current due to Ekman flow, is a function of the wind stress 
and can be estimated by, 


V.| = [| = 1.57 cm sec! = 0.0565 km hour” | (21) 
hf 





O 


where tals the average wind stress, 0.66 dyn cm-*, 0 1s the density of 


Sea water, h 1s average mixed layer depth, 5000 cm and f is the Coriolis 
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parameter, 2wsing, where w is the angular velocity of the earth and ¢ 
is the latitude. 

The velocity calculated due to Ekman flow is insignificant when com- 
pared with the drift rate velocity calculated for R/P FLIP (V=9.3 km hour” !) 
from 28 Jan 1974 to 14 Feb 1974. The calculated Ekman velocity is not suf- 
ficient enough to explain the drift of POLE. Hence, R/P FLIP was by no 
means flowing with the same mixed layer water mass throughout the experi- 
ment. Therefore, the velocity vector V will depend on the drift rate of 
R/P FLIP, V = V.. What Simpson and Paulson {1979} call advection was not 
actual advection for the area but apparent advection due to the drift of 
er FLIP. 

2. Calculation of the Rates of Advection 

The distance R/P FLIP traveled from 0700 hours 28 Jan 1974 to 
0700 hours 14 Feb 1974 (408 hours), was calculated to be 121 km, using 
a rhumbline distance calculation and a position of 35°36'N, 155°05'W to 
34°36N, 155°25'W. Vertical profiles of temperature and salinity were 
taken from R/P FLIP from 0000 hours 30 Jan 1974 through 0700 hours Feb 
1974 (total of 367 hours) which covered a distance of approximately 
109 km. Figures 19 and 20 show the observed changes in heat and salt 
fluxes. The amounts of advection in a coordinate system moving hori- 
zontally with R/P FLIP were then estimated by subtracting the accumu- 
lated surface heat and salt fluxes from the changes in heat and salt 
content, respectively. Figures 21 and 22 show the results of this cal- 
culation. Least square fits were performed on the results, giving a 
Slope of 2.0 cal cm™Chour™ | for heat content gain and a slope of 


1 


0.45 [700 cm hour. for salinity gain. These values can be used to 
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compute the average horizontal temperature and salinity gradients by 
multiplying by the total hours of observation (367 hours) and dividing 
by the average mixed layer depth (approximately 40m) and the distance 
FLIP traveled during the observational period (109 km), 


1 


|VyT| = 0.0017 C km” (22) 


(23) 


Iv4S| = 0.0004 °700 km” 
where }VuT| is the horizontal temperature gradient and [VS | is the hori- 
zontal salinity gradient. 

For comparison purposes, gradients were obtained from a NORPAC 
Atlas [1960] by dividing the atlas temperature change (AT = 0.2 C) and 
the atlas salinity change (AS = 0.05 0700) between FLIP starting and 
final positions (35°39'N, 155°05'W to 34°36'N, 155°25'W) by the total dis- 
tance (121 km) which FLIP traveled, 
[VyTI, = 0.0017°C km”! 
un 


where IVT] a and IVuS| are the horizontal temperature and salinity 


0.0004 9/00 km7! 


gradients that were calculated from the NORPAC Atlas. Although these 
Values agree with the earlier calculation, equations (22) and (23), they 
were only obtained to show agreement in order of magnitude with recorded 
atlas values. These results further illustrate the fact that what 
Simpson and Paulson [1979] call advection was probably the result of the 
Siemnc Of R/P FLIP. 

Figures 19 and 21 show that the temperature field, which is a 
function of both time and depth, has large fluctuations which can only 
be explained by advection, or the drift of FLIP which is a non one- 


dimensional effect. These large fluctuations dominate the heat budget 
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on the time scales shorter than 100 hours, but over a longer time scale, 
surface flux is equally important (one-dimensional effects). Figure 23 
is a plot of observed sea surface temperature and a plot of the model- 
computed mixed layer temperature, with the model-computed temperature 
slightly offset on the temperature scale for comparison. Comparing the 
two plots, a one-dimensional component can be detected in the observed 
sea surface temperature time series but the stronger component due to 
the drift of FLIP is evident. Therefore, this apparent advection (as 
seen by FLIP) of temperature is important on all time scales, but domi- 
nates the short time scale. Figures 20 and 22 show that advection of 
Salinity dominates the salinity budget over both the short and long time 


scales up to at least 400 hours, which was the period of observations. 


D. MIXED LAYER DEPTH 
1. Observed Mixed Layer Depth 

The observed mixed layer depth, shown in Figure 1/7 as used by 
Simpson and Paulson l1979], is the shallowest depth at which the den- 
Sity 1S not more than 0.02 O, units greater than the density at a depth 
of 5m. However, salinity and temperature changes can occur above 
Z = 5m. In such cases, the above definition will result in a misleading 
computed time series for mixed layer depth. A more practical definition 
for observed mixed layer depth, which would be more consistent with the 
precision of the available POLEX observational data, would be the shal- 
lowest depth at which the observed value of on 15002 Oy units greater 
than the observed surface value. For comparison with the Garwood model, 
the observed mixed layer depth (Ao) will be defined according to this 


latter definition. 
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Contours of temperature for the observation period are shown in 
Figures 24 and 25. Figure 24, which is based on a contour interval of 
0.5 C, shows a mixed layer depth as varying from near the surface to ap- 
proximately 90 meters. Figure 25, which has a contour interval of 0.05 C 
after having been smoothed twice,© Shows a mixed layer from near the sur- 
face to approximately 80 meters. This illustrates the importance of 
resolution when determining mixed layer depth from temperature contours. 
A higher resolution frequently results in a shallower mixed layer depth. 
Therefore, a temperature contour interval of 0.5 C is not sufficient 
resolution to accurately determine mixed layer depth from the contour 
plot. Although density is a function of both temperature and salinity, 
during this time series the salinity influence was not significant, so 
that the density was mainly dependent upon temperature. A mixed layer 
depth definition could then be defined as the shallowest depth at which 
the temperature is not more than 0.1 C less than the temperature at the 
Surface which would compare well with the similar definition based on the 
observed density structure, as is evident by comparing Figure 25 with the 


observed mixed layer depth in Figure 26. 


2. Comparison of the Model-Computed Mixed Layer Depth with Observed 
Mixed Layer Depth 


Figure 26 1s a time series of the depth of the model computed 
turbulent boundary layer, -h(t), and of the observed mixed layer depth, 


“ho(t). The value of h should be less than or equal to the mixed layer 


depth as determined from observed 39/dz if the model simulation of the 


é Smoothing was accomplished by giving the center data point a weight 
of 2 and all_of the surrounding data points a weight of 1. 
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vertical density structure changes are accurate. Figure 27 is a plot of 
the same functions, h(t) and h(t) between time = 200 hours and time = 400 
hours. Notice that the model results compare favorably with the observa- 
tions except for a period of about one day starting at hour 285. Figure 
28 is the model computed mixed layer depth, h,(t). The value of -h, (t) 

is taken to be the depth at which the model Of is 0.020, units creater than 
the value of <O4> where < > denotes the vertical average for the turbu- 
lent boundary layer. This definition is equivalent to that for hg{t) ee 
the turbulent boundary layer is homogeneous. Figure 29 shows both h,(t) 
and the observed mixed layer depth, ho it) for the same time period as 
Figure 28. Notice that h)< h for all t. Although hy(t) compares well] 
with h(t), there are times during which hy it) is more comparable to h(t). 
For example, for the period between t = 285 hours and t = 310 hours, 

hy (t) is most similar to h(t). This happens because, although the tur- 
bulent boundary layer may be quite shallow, the associated transient 
thermocline is not sufficiently strong to influence the values for h, (t) 
and hy {t). 

Only during one time period, at about time = 100, is the observed 
mixed layer depth found to be significantly deeper than either the tur- 
bulent boundary layer, h(t), or the mixed layer depth, h,(t), computed 
by the model. This difference between the predicted and observed mixed 
layer depth could be the result of R/P FLIP drifting through a warm core 
eddy or filament of water mass having more tropical characteristics. 
Figures 21 and 22, which are the graphs of the advective components of 
the heat and salt budgets, show a large increase in the magnitudes of 


both the advection of heat and salt at about time = 100 hours. This 
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large increase in advection of both heat and salt is consistent with the 
hypothesis that R/P FLIP drifted across either an eddy or a filament of 
water from the south of the Subtropical front. The Subtropical water 
mass which lies to the south is warmer and more saline than the Eastern 
North Pacific Central water. 
3. Salinity Effects on the Mixed Layer Depth 

An additional model simulation was conducted in order to deter- 

mine what effect precipitation has on the mixed layer depth. A steady 


gm cem™?sec” | was prescribed over a 12- 


rate of precipitation of 5 x 10 
hour period starting at t = 200 hours. At all other times precipitation 
was assumed to be negligible. This event corresponds in magnitude, at 
least, with the rate of precipitation reported by Simpson [1977] for 
5 February 1974. Figure 30 shows a rapid decrease in the surface salin- 
ity with the injection of the precipitation. It drops from about 34.11 2700 
=o 33.99 9700. This drop in surface salinity causes a rapid rise in the 
Surface salinity flux as shown in Figure 31]. Figure 32 shows the mode] 
prediction of the mixed layer depth with and without the precipitation 
event. The associated increase in downward buoyancy flux, Biw' (gy, due 
to the increase in salinity flux (Figure 31) caused the model mixed 
layer depth to deepen much less rapidly than for the case without pre- 
cipitation. The effect would have been even more apparent if the pre- 
cipitation period had not coincided with strong winds. 

Analysis of Figure 29 shows that the turbulent boundary layer 
depth predicted by the model was approximately 20 meters deeper than the 


observed mixed layer depth at t = 210 hours. This time period coincides 


with that reported by Simpson and Paulson [1979 | as the period of most 
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intense precipitation. These results suggest that, if the amount of 
precipitation were accurately known and included in the model boundary 


conditions, there would be an even better correlation between the mode} 


mixed layer depth and the observed mixed layer depth. 
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IIIT. SUMMARY 


This analysis of 15 days of salinity, temperature and density pro- 
files from R/P FLIP as part of the NORPAX POLE experiment and the ac- 
companying model study suggest that, non one-dimensional effects were 
very important for both the heat and salt budgets. Vertical mixing and 
mixed layer depth are controlled primarily by the one-dimensional response 
of the turbulent kinetic energy budget to the local atmospheric forcing. 
It is apparent that the non one-dimensional effects may be largely ex- 
plained by the drift of R/P FLIP. Results from POLE site showed that this 
apparent horizontal temperature advection (as observed by the moving plat- 
form) was important on all time scales and dominated the smal] time scale. 
Nevertheless, the one-dimensional model computation of mixed layer depth 
compared favorably with the observed mixed layer depth. It is interest- 
ing that even though this apparent advection usually dominated the salt 
budget over both the long and short time scales, the one-dimensional mixed 
layer model response to the surface salinity flux can be significant. 
During the single period of strong precipitation, the mixed layer depth 
was shallower (on the order of tens of meters) than it would have been 
without precipitation. This is consistent with the results of Miller 
[1976] , showing that salinity flux can at least occasionally play an 
important part in mixed layer dynamics in the subtropics. Further 
Studies would have to be conducted with observational data that included 
the time, duration, and amount of precipitation before a firm conclusion 


can be made on the longer-term importance of salinity effects on the 
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mixed layer depth. Qualitatively, model performance was improved when 
the influence of precipitation was taken into account. 

Simpson and Paulson [1979] concluded that the one-dimensional models 
failed to predict the mixed layer depth and temperature observed during 
POLE. Why did the models of Pollard et al [1973] and Niiler [1975] , which 
were used by Simpson and Paulson, fail to predict the observed mixed layer 
depth? One part of the explanation for the model's apparent failure could 
be the definition of mixed layer depth used by Simpson and Paulson. An 
accurate determination of the mixed layer depth from a vertical profile 
can only be made by considering the entire profile, an erroneous analysis 
of the mixed layer depth is probable whenever there is a transient 
pycnocline above z = 5 meters. Therefore, contrary to the results of 
Simpson and Paulson [i979] , it has been shown that the one-dimensional 
turbulence closure model used here does seem to explain most of the ob- 
served mixed layer depth changes and that the mixed layer depth is more 


influenced by the local surface fluxes than by horizontal advection. 


34 





7 ~~ 
PREVAILING WESTERLIE 
a <<” ' 
D> > OR7P FLIPZ 
, oe rN 
Ue LTT: VM 
REGION OF SUBTROPICAL FRONT 


oI LiL Li ase VU Re < 


25° aa 
i _ REGION 


BD” 





Figure 1. The location of the R/P FLIP during the POLE 
experiment (35°N, 155°W) in relation to general 
oceanic features. (Simpson and Paulson 1977) 
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Figure 2. 


Idealized density and mean velocity profiles of the ocean 
mixed layer (----). Mixed layer depth is (h), (6) is the 
thickness of the interface or entrainment zone. Region I 
is the fully turbulent mixed layer depth. Region II is the 
slightly stable, intermittently-turbulent entrainment zone. 
Region III is the stable underlying watermass having negli- 
gible vertical fluxes in comparison to those of Region I. 
(Garwood 1977) 
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Figure 3. 
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Typical POLEX data temperature, salinity, and sigma-t pro- 
files taken from Simpson and Paulson [1979]. Added to the 
profile is Region I the fully-turbulent mixed layer of 
depth h. Region II the slightly stable, intermittently- 
turbulent entrainment zone of thickness 6. Region III is 
the stable underlying water mass having negligible fluxes 
in comparison to those of Region I. 
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Figure 4. 
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Typical POLEX data temperature, salinity, and sigma-t pro- 
files taken from Simpson and Paulson [1979]. Added to the 
profile is Region I the fully-turbulent mixed layer of 
depth h. Region II the slightly stable, intermittently- 
turbulent entrainment zone of thickness 6. Region III is 
the stable underlying water mass having negligible fluxes 
in comparison to those of Region If. 
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Figure 17. Interpolated hourly values of mass and buoyancy fluxes, 
surface salinity, sea surface temperature, surface 
Sigma-t, and the various components of the heat flux 
involved in the determination of the buoyancy flux. 
(Simpson and Paulson 1977) 
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Figure 18. Interpolated hourly values of the various components 
of the surface momentum and heat balance are shown. 
(Simpson and Paulson 1977) 
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Figure 19. 
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The observed changes in heat content in the upper 
40 meters (solid line). Model computed accumulated 
surface heat flux (dotted line). The origin of the 
abscissa is 0000 hour (local) 30 January 74. 
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The observed changes in salt content in upper 40 
meters (solid line). Model-computed accumulated 
Surface salinity flux (S'wi), dashed line, which 
is negligible in comparison with the observed 
changes in salt content. The origin of the 
abscissa is 0000 hour (local) 30 January 74. 
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Advection of heat (solid line) determined by © 
Subtracting the accumulated heat flux from the 
change in heat content (Figure 21). Least 
squgres rit to data resulted in a slope 2.0 cal 
cm ~ hour 
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Advection of salt (solid line) determined by 

subtracting the accumulated surface salinity 

flux from the changes in salt content (Figure 
20). Least square fit to data resulted ina 

slope 0.45 0/00 cm hour-!. 
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Figure 24. Isotherm distribution for the period 30 Jan 74 through 
14 Feb 74. The contour interval is 0.5°C. 
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(0) is the observed mixed layer depth, 


The solid line is the model computed turbulent boundary layer 


The origin of the abscissa is 1000 GMT 28 January 1974. 
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